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The dilational viscosity e' of the Langmuir monolayer is considered in a the- 
oretical model taking into account an orientational effect of the dilational 
wave on the surface molecules. The surface tension treated in the frame- 
work of this model is supposed to be dependent on the degree of molecular 
ordering in the monolayer plane. This orientational order is described by 
the surface order parameter Q and the orientational part of the free energy 
which is given by Landau's expansion in powers of Q. The magnitude of sur- 
face viscosity determined by the surface tension derivative dj/dQ is found 
to be in good accord with the experimentally observed e'. The sign of e' 
is positive which indicates that increased ordering in the monolayer plane 
increases the surface tension. 

PACS: 61.25.Mv; 68.10.-m; 68.10.Cr 

Strong damping of capillary waves on the monolayer-water interface has 
been well known for many years. However the physical nature of the two sur- 
face viscosities which are included into a dispersion relation |IJ ^] describing 
the propagation of surface waves still remains unclear. It has been suggested 
that the monolayer possesses up to three surface viscosities: dilational, shear 
in the interfacial plane and shear normal to the surface plane [|], |3|]. Only 
two of them: dilational e' and transverse shear 7' do couple to surface modes 
which scatter light These surface viscosities were introduced as the re- 
sponse functions which appear after expansion of the complex moduli 7, e, 
to linear order in u> |jj . Hence the dilational modulus is usually written as a 
superposition of the dilational surface viscosity e' and the surface elasticity 
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e = e + iuie' (1) 

The surface viscosity e 1 so introduced is a pure phenomenological property 
which was suggested as the surface excess quantity || [3| . 

Only for two-component systems of soluble surfactants is there a strong 
theoretical basis for the dilational surface viscosity gj. In such a system 
the surface viscoelasticity appears as a result of the competition between 
adsorption and diffusion in the concentrational boundary layer beneath the 
surface. The main result of this formalism is that e' is not an independent 
property but is a function of some relaxation time r which depends on the 
diffusion coefficient D and the slope of the adsorption isotherm ^ ]5|, |6| : 

e > = e «L ( 2 ) 

where r = ^-^D Juj, (C is the bulk concentration and T s the surface con- 
centration). However for an insoluble (spread) monolayer on a liquid sub- 
strate there is no adequate theoretical description for the dilational surface 
viscosity. At the same time the experimental findings of negative e' in some 
liquid systems J?], |8| indicate that some important physical processes are not 
accounted for in the theoretical formalism of |l|, |(| . 

Our previous papers concerned the dilational viscosity of various liquid 
surfaces ^,10]. It was demonstrated that the surface dilational modulus 
can be described by relaxation of fluctuations of temperature and surface 
charge density. Unfortunately the dilational viscosity calculated from those 
models appears to be a few orders of magnitude smaller than that observed 
in experiments on different liquid systems [0, || . The present paper demon- 
strates that the orientational effect of the surface wave on long molecules can 
be the reason for the dilational viscosity of the liquid-monolayer interface. 

The orientational effect of a shear flow on long molecules in the liquid 
bulk has been known for many years [11]. However such an orientation of 
molecules in the surface layer has not yet been reported. We concentrate 
upon a monolayer being in a liquid-expanded (LE) state where the surface 
area per molecule is relatively large according to the surface isotherm |l2]. 
Due to this fact and strong electrostatic interaction between the molecular 
head groups and water the long axes of molecules are not oriented along the 
surface normal but comprise some angle to it. 
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Unfortunately LE state is not too much studied experimentally com- 
pared to a condensed state. Therefore we will use the results of theoretical 
predictions on order parameters in LE state. One of them develops a 
concept of orientational order in a monolayer comprising hard rods grafted 
on the surface. The main result is that molecular alignment increases contin- 
uously with increasing surface concentration. Since the molecular azimuths 
are not arbitrary but are defined by the tilt angle, the molecular alignment 
in the surface plane is induced by collective tilt. A numerical analysis of 
flexible molecules [O] shows that in-plane ordering in a liquid phase can 
be even enhanced by chain length. All these papers presume that surface 
fluctuations destroy the in-plane long range order. We will give a highlight 
to the effect of dilational waves on the surface order. As it will be shown 
below a very small magnitude of the order parameter is sufficient for our 
model. 



As mentioned in [13|, a monolayer formed by amphiphile diblock copoly- 
mer on water surface can be a realization of the model describing alignment 
of mobile rods. Another experimental situation relevant to the present con- 
text is monolayer of a nematogen on water studied in fl5| . At low surface 
concentrations nematic molecules lie on water but the molecules interact 
with each other as in a 3-D nematic. 

In the present context we suggest that projections of molecules display 
a quasi long-range ordering in the surface plane. In the case of a high tilt 
angle the short axis of a molecule nearly coincide with its surface projection. 
Let us consider the surface footprint of a tilted molecule (see Fig. 1). The 
long axis of the molecular "cross-section" of a length L makes an angle 4> 
with wavevector q (supposed to be along the x-axis) of the dilational wave. 
This molecule experiences an orientation torque induced by the longitudinal 



gradient du x /dx of velocity in the dilational wave as it was suggested in 16]. 
If the center of the molecule moves with the liquid then the two ends will 
move with different relative velocities u\ and U2- 

u\ — U2 = 5u = —r~L cos <fi (3) 

ax 

The transverse (to the molecule axis) component of the relative velocity 
gives rise to rotation of the molecule with an angular velocity f2 = <fi 

ft = _^£ sin2 ^ (4) 
ax 
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where negative sign means that the molecule rotates towards to decrease in 

The angular velocity of the molecular rotation is proportional to the 
gradient of surface velocity whose magnitude is small and oscillating: u ~ 
u'exp(i(u;i — qx)) . The effect is superposed with a strong anisotropy of 
nematic in-plane ordering supposed in the present text. The order parameter 



tensor can be introduced in the pure 2-D case, analogously to 3-D case 17], 
as : Qij = Q(niTij — \5ij) where i=l,2, nj is the nematic director lying in the 
surface plane. Due to the surface wave the order parameter tensor comprises 
two parts: steady Q°- and an oscillating one 

In order to find out the equation describing the behavior of Q' we use 
the well-known phenomenological equation [|17]] : 

where F is the free energy density, A is a proportionality coefficient, Sik is 
the tensor of viscous stresses and \x is some coefficient having the dimension 
of viscosity. Usually the free energy is expanded in powers of Q and the 
expansion contains the terms proportional to the second, third and fourth 
power of Q (see for example (T^j). The main feature of the 2-D case is an 
absence of a cubic term in the energy expansion due to the symmetry in the 
surface plane. Thus the free energy density is written as: 

F = AQijQij + BQijQijQkiQki (6) 

where A is the temperature dependent coefficient changing sign at the point 
of an isotropic-nematic transition and B> 0, and the summation is taken 
over repeated indexes. In principle, it might be possible to calculate A 
and B in the 2-D case as it is done in the bulk nematic [pC[| . However, an 
anisotropic interaction potential is essentially unknown to the present study. 
Therefore the coefficients A and B will be kept without explicit expressions. 
Since we have only one non-zero component of the stress tensor :sn = 



f Tte ( see 0) E( l( 6 ) is simplified 



dt ~ 3 dx rdQu [ ' 

where f is dimensionless energy density / = F/Fq/Vq and r = [xFq/Vq 
is some relaxation time. The derivative dF/dQij can be found from the 
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substitution of the averaged plus the disturbed component of the order pa- 
rameter tensor into EqQ. However we simplify the situation suggesting 
that the molecular axes are oriented in the surface plane either parallel to 
the wavevector q of the dilational wave or perpendicular to it. Then the 
disturbed component of the order parameter tensor has a diagonal form, i.e. 
only two non-zero components: Q' n and Q' 22 whereas Q22 = "Qn = 0! ■ Thus 
the free energy derivative is 
Bf 8 

nk = aQ' {A(Q " + Q ' ? + m ' , + Q ' ),) (8) 

where Qo is the steady part of the order parameter and Q' is the magnitude 
of an oscillating part. Then Eq(^) can be re-written to the first order in Q' 

^ = i^Q< = - l -Q'{2A + UBQl) + (9) 

A is a phenomenological coefficient such that the term A^r = (Q) expresses 
the mean torque exerted on the order parameter by the dilational wave. 
Following Ref |jl9| this is plausible to suppose the coefficient A is the magni- 
tude of the mean angular velocity averaged over all molecular orientations: 
The averaging gives A = 1 if the angular velocity in the form of Eq(||) is 
used. Due to its phenomenological origin, A is not inevitably equal to 1 but 
may depend of molecular shape, or diffusion coefficient, therefore only an 
approximation A ~1 will be used. Hence Q' is 

Q' = -^ - (10) 

^ 3 dx iujT + 2A + 6BQl v ' 

The amplitude of Q' is to dependent if the typical wave frequency and 
relaxation time are such that: u ~ r. We remind that only the real part 
of the complex Q' is physically meaningful; a phase of Q' corresponds to 
some angle between the wavevector q and the averaged orientation of the 
director. The effect of this phase shift is an additional wave damping on the 
surface possessing some anisotropy. Our situation strongly resembles the 
propagation of longitudinal acoustic waves through an anisotropic liquid 



which was considered by Frenkel in [11|. He demonstrated that the com- 
ponents of the order tensor are proportional to the longitudinal gradient of 
velocity in a sound wave. If the wave frequency is comparable to the charac- 
teristic relaxation time then the order tensor and the velocity gradient are 
connected via the usual Maxwell formula describing the viscoelastic relax- 
ation of molecules subjected to oscillating disturbances. Such viscoelasticity 
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points to the existence of a phase shift between a disturbance and the media 
response 
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Our goal is to find the surface elastic modulus which can be written as 

§: 

d~f <97 dQ 

dp s dQ dp s 

where 7 is the surface tension and p s is the surface concentration. The 
linkage between disturbances of the surface concentration and the surface 
velocity is found from the conservation of matter at the surface H || : 



dp s _ n d 2 p s _ du x 

o, — Da 9 PsO 
Ot OTT OX 



Ds^-Pso^ (12) 



where p S Q is an equilibrium surface concentration. We are looking for the 
surface concentration in the form of a wave: p' s ~ g^M-?^) _ Substitution in 
Eq(12) yields: 

PsQ ~dx = p ' s ^~ iuJ + Ds 1 2 ) ~ ~ iuJ P' ( 13 ) 
where we supposed the surface diffusion coefficient D s to be nearly equal to 
that of the bulk (D ~ 10 -5 cm 2 /s) and typical wave frequencies are such 
that uj » D s q 2 gg. It is easy to see that the second derivative in Eq(|Tl|) 

= JL ^ (14) 

d Ps Z Ps0 (iuT + 2A + l2BQl) y ' 

where a negative sign means that Q is a decreasing function of the surface 
concentration. This is a typical signature of surface fluctuations; they tend 
to destroy a long range nematic order whose order parameter usually grows 
with p s in the absence of fluctuations |l3| . 

Some papers and references herein) postulate the existence of a 

nematic ordering induced on the surface of a liquid crystal being in an 
isotropic state having Q=0. This surface ordering treated in [^] is due to 
the contact with a solid wall inducing a molecular alignment rapidly decaying 



into the liquid bulk. In [21] the free surface energy has been written using 
the Landau expansion over the in-plane order parameter. Note that liquid 
crystals have relatively large penetration depth of the surface ordering while 
the ordering in spread monolayers is strictly confined to the first surface 
layer. Therefore the difference from liquid crystals is an absence of the 
Frank elastic modulus describing the reaction of the bulk on the surface 
ordering. 
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We establish now an analytical dependence of the surface tension on 
the surface order parameter. For the Langmuir monolayer only the surface 
has some ordering-the liquid bulk is essentially isotropic. Therefore the 
orientational part of the surface free energy is given by the integration of 
F over the liquid depth: F s = J °° F(Q)5{z)dz = £F(Q S ), where S(z) is the 
delta function and Q s is the magnitude of the order parameter on the surface 
and £ is some correlation length. Thus the orientational part of surface free 
energy density is written as: 

F s = ik B TN APb /M m {AQ 2 s + BQ A S ) = k B TN A -^F(Q), (15) 

lVl m 

where M m is the molecular weight and the relation p s = connects the 
number of particles in the surface and the bulk. It should be borne in 
mind that only a minimum of the function F(Q) corresponds to physically 
stable state |l7| , p2| and minimization with respect to Q gives: F rn i n =A 
Qo/2, where Qq=-A/2B [^]. It is apparent that in the 2-D case we have a 
continuous phase transition instead of the first order transition in the bulk 
liquid crystal. 

The derivative djo/dQ can be found from a general thermodynamic 
equation relating the surface tension and the surface excess density of mass 
r s @: 

dj = -T s dfi s (16) 

where fi s is the surface chemical potential. Note that the entropy term - 
S s dT is omitted due to isothermal character of our effects. The chemical 



potential is fj, s = (dF s / 8N)a, by the definition [22] and T s = p s for any 
insoluble film. The derivative d~f/dQ is 

dQ = ~ Ps dQ = - pskBT -dQ- (17) 

The variations in the surface mass density upon a phase transition are neg- 
ligible. 

As we mentioned above only a minimum in the free energy is realized 
in the equilibrium. One might think that dF/dQ is zero as it should be in 
the minimum of a function. However this is not so due to deviations from 
the equilibrium induced by surface wave. To calculate dF/dQ one should 
take into account the changes in surface area (and surface pressure) due to 
dilational wave. This means that the position of a minimum in the free 
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energy and the corresponding order parameter Qo also would change upon 
surface compression (rarefaction) due to a strong pressure dependence of A. 
Therefore we can take the Qo derivative, having in mind that the position 
of minimum shifts with the change in A: dF/dQo = —AQq. This derivative 
is negative (Qo > 0, A< in nematic phase) and the minimum of the free 
energy is lowered upon an increase in Qo- Hence the Q derivative of the 
tension is positive: 

d 1 /8Q = Ps /M m N A k B T \A\Q (18) 

It turns out that the derivative d^/dQ is positive for the whole Q- 
range and the surface viscosity is also positive according to Eq(|iT|) where 
dQ/dp s <0 . Searching additional physical arguments in favor of an increas- 
ing 7(Q) dependence one can return to Eq(|ll"D introducing surface elasticity 
in the framework of orientational model. We may perturb a film by a small 
horizontal displacement which rarefies locally the monolayer and decreases 
the surface concentration. At the same time the surface expansion would 
increase the surface order parameter due to the orientational effect of the 
dilational wave discussed above. Negative sign in the definition of e re- 
flects the feed-back reaction of a perturbed surface element. For instance, 
one may expect the surface tension to decrease on expansion of the surface 
which is accompanied by some increase in Q. In this case a restoring force 
would act outward i.e. from the surface region with lower 7 to those regions 
having higher tension. Thus a decreasing 7(Q) dependence would lead to a 
semi- infinitive expansion and the film instability. On the contrast, having 
supposed an increasing 7(Q) dependence the surrounding film would tend 
to compress the surface element which was expanded initially. This ensures 
that increasing 7(Q) provides the monolayer stability setting up d'y/dQ >0 
and positive e'. 

The modulus of dilational viscosity given by e' = Q[e]/u, found from 
Eqs. ©, (0), © 

// s _ k B TN APs 1 A I Qq it 

H) 'u\M m it + 2/uj(A + 6BQq) 

is plotted in Fig 2 for different time constants t using a typical value of the 
surface concentration p s ~ 1 mg/m 2 . This plot ensures that the dilational 
viscosity exhibits a strong uj dependence for some r > 10 -5 s while for 
shorter time constants e' is nearly independent of uj. 
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Surface viscoelastic properties of glycerol mono-oleate monolayers on wa- 
ter were studied by surface light scattering in p3fl . Surface properties were 
deduced from a many parametrical fit to the measured autocorrelation func- 
tions. Although the surface viscosity was introduced in the framework of 
a phenomenological model (see above) e' ~ 10~4 mN s/m was found. The 
dilational viscosity inferred from the fit was sensitive to the monolayer com- 
pression: the increase in ef about an order of magnitude upon the monolayer 
compression have been reported. 

Our analysis shows that the dilational surface viscosity appears in the 
form of oscillations (with some relaxation time) of the order parameter which 
take place due to the passage of the dilational wave. This surface anisotropy 
provides a phase shift between the velocity gradient and the in-plane direc- 
tor. The surface dilational viscosity introduced in our model is positive and 
drops to zero at the point where A>0 and Q=0 i.e a point of the surface 
phase transition of a second order. The main assumption of the present 
paper is the existence of nematic ordering in the liquid-expanded state. 
However the present model of surface viscosity can be modified for hexatic 
ordering in Langmuir monolayers. 

Fruitful discussions were shared with E.I. Kats. 
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Captions 



Fig.l Orientational effect of the velocity field in the surface wave on the cross 
section of a molecule in Langmuir monolayer, q is the wavevector of 
the dilational wave. The dashed lines show the lines of equal phase of 
velocity in longitudinal wave. Surface velocity remains constant along 
the lines but varies perpendicularly to them giving rise to the gradient 
du^/dx on the surface. 

Fig.2 Surface dilational viscosity e'(u) calculated just near the surface phase 
transition for different relaxations times r. r are labeled near each 
curve in s. The following values of the parameters are used: Q=0.1, 
| A |=0.4, B=l 
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